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ABSTRACT 

We present a simplified and fast method for simulating minor mergers between galaxy 
clusters. Instead of following the evolution of the dark matter halos directly by the 
N-body method, we employ a rigid potential approximation for both clusters. The 
simulations are run in the rest frame of the more massive cluster and account for 
the resulting inertial accelerations in an optimised way. We test the reliability of this 
method for studies of minor merger induced gas sloshing by performing a one-to-one 
comparison between our simulations and hydro+N-body ones. We find that the rigid 
potential approximation reproduces the sloshing-related features well except for two 
artefacts: the temperature just outside the cold fronts is slightly over-predicted, and 
the outward motion of the cold fronts is delayed by typically 200 Myr. We discuss 
reasons for both artefacts. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: A2029 X-rays: 
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1 INTRODUCTION 

During the last decade, high resolution X-ray observations 
have found a wealth of structure in the intra- cluster medium 
(I CM) of galaxy clusters, among them cold fronts (see re- 
view by Markevitch k Vikhlinin 2007). These structures re- 
veal themselves as sharp discontinuities in X-ray brightness 
accompanied by a jump in temperature, where the brighter 
side is the cooler one. One variety of cold fronts (CFs) was 
soon understood to be the contact discontinuity between the 
gaseous atmospheres of merging clusters (e.g. A2142: Marke- 
vitch et ah 2000; A3667: Vikhlinin et al. 2001; and the bullet 
cluster IE 0657-56: Markevitch et al. 2002). 

A second class of CFs was found to form arcs around 
the cool cores of apparently relaxed clusters (e.g. RX 
J 1720.1-1-2638: Mazzotta et al. 2001; Mazzotta &; Giacin- 
tucci 2008; Owers et al. 2009; MS1455.0+2232: Mazzotta k 
Giacintucci 2008; Owers et al. 2009; 2A0335+096 Mazzotta 
et al. 20C3; Sanders et al. 2009; A2029: Clarke et al. 2004; 
Million k Allen 2009; A1795: Markevitch et al. 2001; Bour- 
din k Mazzotta 2008; Perseus: Churazov et al. 2003; Sanders 
et al. 20C5; A496: Dupke et al. 2007; Ghizzardi et al. 2010; 
Virgo: Simionescu et al. 2010; Centaurus: Fabian et al. 2005; 
Sanders k Fabian 2006). Markevitch et al. (2001) suggested 
that this variety of CFs forms due to sloshing of the cool 
central gas within the central cluster potential, where the 
sloshing is initially triggered by a minor merger event. Using 
hydro+N-body simulations, Ascasibar k Markevitch (2006) 
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(AM06 afterwards) have shown that the sloshing scenario 
reproduces the morphology of observed CFs. 

In recent years, more and more sloshing CFs have been 
reported (e.g. Owers et al. 2009; Ghizzardi et al. 2010). 
In principle, the properties of the CFs contain informa- 
tion about the merger history of the clusters, as has been 
shown for the Virgo cluster (Roediger et al. 2011, Rlla 
hereafter), A496 (Roediger et al. 2011, Rllb hereafter), and 
RXJ1347.5-1145 (Johnson et al. 2011). However, disentan- 
gling the merger history of each cluster from the CF proper- 
ties requires a set of dedicated simulations for each cluster. 
Doing this with full hydro+N-body simulations is computa- 
tionally expensive. The same is true if the influence of more 
time-consuming physics like viscosity (ZuHone et al. 2010), 
magnetic fields (ZuHone et al. 2011) or thermal conduction 
(ZuHone et al. 2011, in prep.) on sloshing CFs is studied. A 
reasonable simplification that speeds up the simulations con- 
siderably would be very useful. The most expensive part of 
such simulations is the seif- gravity of the gas and dark mat- 
ter (DM) particles. In most cases, the self-gravity of the gas 
can be neglected, because the Jeans length is about 1 Mpc. 
The simulations speed up substantially when the DM ha- 
los of the main cluster and the subcluster are approximated 
as rigid potentials (RPs). However, AM06 have shown that 
even in a minor merger the central part of the main cluster 
moves significantly w.r.t. the overall cluster potential, thus 
this effect cannot be neglected. ZuHone et al. (2010) (Z10 
hereafter) suggested that the rigid potential approximation 
can be used if additionally a point-mass-like approximation 
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of the motion of both clusters is taken into account, includ- 
ing the effects of inertial acceleration. 

Here we improve this rigid potential approximation 
and demonstrate the reliability of the simplified simulations 
for the application to gas sloshing by comparing them to 
hydro+N-body simulations. This rigid potential approxima- 
tion has already been applied successfully in sloshing simu- 
lations for the Virgo cluster (Rlla) and for A496 (Rllb). 


2 METHOD 


2.1 Test setup 


We consider the following scenario: The I CM in a mas- 
sive spherical galaxy cluster (the main cluster) is initially 
in hydrostatic equilibrium. A spherical, less massive, gas 
free galaxy cluster (the subcluster) passes through the main 
cluster. The gravitational impact of the subcluster initiates 
sloshing of the ICM in the main cluster core and subsequent 
cold front formation (see AM06 for a detailed description of 
the dynamics). 

Hydro-fN-body simulations of this scenario have been 
performed by AM06 and Z10. Our aim here is to investi- 
gate to what extent simulations with a rigid potential ap- 
proximation (RPA) for the DM content of, both, the main 
cluster and the subcluster can reproduce the resulting CF 
structures in terms of morphology, orientation, size, tem- 
perature and density distribution. As the reference, we use 
the hydro-|-N-body simulations of Z10. Thus, we follow their 
lead and tailor our initial models to match theirs. The main 
cluster model is based on a Hernquist potential (Hemquist 
1990) with a scale radius, a. The temperature profile is de- 
scribed by the phenomenological function 


tm - T ° c + r / g ° 
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where To is a measure for the overall cluster temperature, c 
describes the depth of the central density drop, and a c char- 
acterises the radius of this drop (see also AM06). The corre- 
sponding density profile resulting from hydrostatic equilib- 
rium is (AM06) 
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In our setup, we initialise the main cluster with these den- 
sity and temperature profiles, choosing the parameters such 
that they fit the corresponding hydro+N-body simulation 
(see Table 1). From these profiles, we derive the gravita- 
tional potential of the main cluster assuming hydrostatic 
equilibrium. 

As in Z10, the subcluster is initialised as a pure DM 
structure, where the DM mass distribution is described by 
a Hemquist profile (Hernquist 1990). 


2.2 The rigid potential approximation 

Our simulations are run in the rest frame of the main cluster. 
The gas dynamics is described by the hydrodynamical equa- 
tions. Additionally, the ICM is subject to the gravitational 
acceleration due to the main cluster and the subcluster. 


2.2.1 Orbit of the subcluster 

We assume the orbit of the subcluster to be the orbit of a test 
particle free-falling through the main cluster. In the course 
of the simulation, the potential of the subcluster is shifted 
through the main cluster along this orbit. This approach 
does not include dynamical friction, which will slow down 
the subcluster after pericentre passage (see Sect. 3.3 for a 
comparison). However, the sloshing is triggered mainly dur- 
ing the pericentre passage, and thus our results do not suffer 
from this discrepancy. Given that our test particle orbit is 
bound to predict increasingly wrong subcluster positions af- 
ter the first pericentre passage, we stop our simulations well 
before a second pericentre passage. 

We construct orbits that are comparable to the ones 
in the hydro+N-body simulations in orientation, pericentre 
distance, and velocity history prior to pericentre passage. We 
start our simulations 1 Gyr prior to the pericentre passage 
of the subcluster, which we set to occur at t = 0. We reran 
our fiducial simulation with an initialisation time of 0.5 Gyr 
prior to pericentre passage and found that our results are 
not sensitive to the choice of the initialisation time. 

The subcluster orbit is in the zy-plane of the computa- 
tional grid. For the sake of a short notation, we identify the 
+y-direction as ’’north” (N), the — ^-direction as 55 south” 
(S), the +x-direction as ’’west” (W), and the — ^-direction 
as ’’east” (E). The subcluster will start W of the main clus- 
ter core, has its closest approach to the main cluster core in 
the NE and moves away towards the SE. 


2.2.2 Basic inertial frame correction 

As the rest frame of the main cluster is not an inertial frame, 
the ICM in this frame is subject to a pseudo-acceleration due 
to the attraction of the main cluster core towards the ap- 
proaching subcluster. Z10 used the most simple approxima^ 
tion to account for this: they assumed that the main cluster 
responds to the gravity of the subcluster as a whole, like a 
rigid body. Consequently, they calculate the inertial accel- 
eration felt by the main cluster centre due to the subcluster 
and add this pseudo-acceleration to all of the ICM. 


2.2.3 Improved inertial frame correction 


This basic inertial frame correction is a reasonable approx- 
imation for the central region of the main cluster, but it is 
wrong for the outer parts of the cluster. It will lead to unreal- 
istic flows in the outer parts of the main cluster. For smaller 
pericentre distances, during pericentre passage of the sub- 
cluster the inertial acceleration can be large and even pro- 
duce supersonic motions in the cluster outskirts. These un- 
realistic flows can influence the resulting CFs at later stages. 
Hence, we propose to apply the pseudo-acceleration only to 
the central region of the main cluster inside a characteristic 
radius, i?damp, and dampen it outside this radius exponen- 
tially over a length scale, Ldamp. Thus, instead of adding 
the same inertial frame acceleration at every position in the 
cluster, we multiply it with a radius-dependent function, 


W(r) = 


{ 


i 


exp(- 


r -^damp 

^damp 


if r < Rdamp 
else 


( 3 ) 
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Table 1. Summary of model parameters. See Sect. 2.1 for details. 


merger characteristics 


mass ratio 

20 

5 

2 

impact parameter 

200 kpc 

500 kpc 

500 kpc 

pericentre distance 

62.5 kpc 

150 kpc 

135 kpc 

main cluster 

/°o/(10 -25 gem -3 ) 

3.75 

2.9 

2.76 

a/ kpc 

615 

623 

615 

n 

5 

5 

5 

2o/(10 7 K) 

15.8 

13.8 

11.27 

fle 

61.5 

62.3 

61.5 

c 

0.1585 

0.169 

0.1557 

subcluster 

mass (10 14 Mq) 

0.714 

2.5 

5 

scale radius (kpc) 

220 

350 

416 

damping 

none, 

none, 

none, 

(Rdamp/ kpc, 
Ldamp/ kpc) 

(500, 300) 

(800, 300), 
(800, 100), 
(500, 300) 

(950, 400) 


where r is the distance to the main cluster centre. 

The choice of Rdamp is motivated by the sphere of in- 
fluence of the subcluster: If the subcluster passes the main 
cluster centre at a pericentre distance larger than its own 
size, it attracts the main cluster core only slightly. If the 
subclustei passes the main cluster core at a small distance, it 
can attract at most a region comparable to its own size, but 
not beyond that. Hence, Rdamp should be comparable to the 
size of the subcluster, i.e., about twice its characteristic scale 
length. We test several combinations of (Rdamp, Ldamp). 

2.3 Code 

All simulations are run with the FLASH code (version 3.2, 
Dubey et al. 2009). FLASH is a modular block-structured 
AMR code, parallelised using the Message Passing Interface 
(MPI) library. It solves the Riemann problem on a Carte- 
sian grid using the Piecewise-Parabolic Method (PPM). The 
simulations are performed in 3D and all boundaries are re- 
flecting. We use a simulation grid of size 4 x 3 x 3Mpc 3 , 
which is large enough to prevent reflected waves reaching 
the central region of interest during our simulation time. 
We resolve the inner 50 kpc with 2 kpc, the inner 130 kpc 
with 4 kpc, the inner 260 kpc with 8 kpc and enforce de- 
creasing resolution with increasing radius from the cluster 
centre to optimise computational costs. We have performed 
resolution tests here and also in Rlla,b and found our re- 
sults to be independent of resolution. 


3 RESULTS 

We present simulations of three merger scenarios which 
cover a range in mass ratios from 2 to 20. The complete 
model parameters for each run are listed in Table 1. 

In the following subsections we give a detailed compar- 
ison of the RP simulations to the hydro+N-body ones. In 
order to aid the reader in smoothly following our analysis, 


we state our main finding already here: The RPA is able to 
reproduce the hydro+N-body simulations very well except 
for two systematic differences, of which one can be fully cor- 
rected for, and the other partially: 

• After the onset of sloshing, the RP simulations lag be- 
hind the hydro+N-body ones by 200 to 250 Myr, depending 
on the cluster mass ratio. This lag can be corrected for by 
delaying the hydro+N-body results by the appropriate lag 
in the comparison. 1 

• The RPA generally over-estimates the temperature just 
outside the CFs. This disagreement can be attenuated by the 
damping discussed above, but not avoided completely. 

3.1 Fiducial case: mass ratio 5, moderate damping 

Our fiducial case is the intermediate merger with a mass 
ratio of 5 between the clusters. We simulate this case 
with four different settings for the inertial frame correc- 
tion as listed in Table 1: one without large-scale damp- 
ing, two moderate damping settings and a strong one. 
The best results are achieved for moderate damping with 
(■Rdamp/ kpc, Ldamp/kpc) = (800, 300). This case will be de- 
scribed first. The effect of the damping will be summarised 
in Sect. 3.2. 

3.1.1 Morphology and orientation of cool spiral 

With both simulation methods, the gas sloshing evolves in 
a very similar manner. We demonstrate this in Fig. ' 1 by 
showing snapshots of the ICM temperature in the orbital 
plane. The first and third row are for the hydro+N-body 
run, the remaining two rows for our fiducial RP simulation. 

The top two rows focus on the onset of sloshing. We 
see the subcluster pass the cluster centre (timesteps 0 and 
1 Gyr) from the NW over NE towards the SE. At 0.5 Gyr, 
sloshing has just set in, and an arc-like CF towards the S is 
accompanied by a cool fan towards N. Both are surrounded 
by hotter ICM. These properties are alike in both methods. 
In the RPA, the cool fan takes a more spiral-like appearance 
compared to the hydro+N-body case. Also the temperature 
distribution S of the southern CF differs between both meth- 
ods. 

The bottom six panels of Fig. 1 display the further evo- 
lution of the gas sloshing. As mentioned above, we find that 
the RP simulations lag behind the hydro+N-body ones by 
250 Myr for this mass ratio. Further details of this lag will 
be discussed in Sect. 3. 1.2. 2. Hence, here we plot the results 
from the hydro+N-body simulations with a delay of 250 Myr 
in order to compare corresponding timesteps. 

In the intermediate phase (0.7 to 1.5 Gyr) the cool spi- 
ral typical for sloshing forms. The major CF is found in the 
SW, and a secondary CF evolves towards the NE. At its 
outside, the cool spiral is surrounded by a hot horse-shoe 
shaped region, which tends to be slightly too hot in the 

1 In a correct manner, the results of the RPA should be brought 
forward instead of delaying the more accurate hydro+N-body 
ones. However, there are several RP runs for each merger case 
and only one hydro+N-body run. For the sake of simplicity and 
clarity we decide to apply the time shift to the hydro+N-body 
run and trust the reader to remember this footnote. 
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Figure 1 . Comparison of hydro+Nbody simulations (first and third row) and rigid potential simulations (second and fourth row) for 
the merger with, mass ratio 5. The panels show the ICM temperature in the orbital plane (see colour scale). The panel size is 1 Mpc. The 
timestep is noted above each panel. The hydro+Nbody results are from Z10, replotted in our colour scale. The rigid potential simulation 
shown here uses the damping setting (Rd &mp / kpc, Tdamp/ kpc) = (800, 300). The upper six panels show the onset of sloshing, the bottom 
six the evolution of the cold fronts and cold spiral structure. Our analysis shows that after the onset the rigid potential simulations lag 
behind the hydro+Nbody one by 250 Myr (see Sect. 3. 1.2. 2). Hence, in the bottom six panels we plot the results from the hydro+Nbody 
simulations with a delay of 250 MyT. 
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RPA. This evolution is the same in all cases, and the results 
of both methods agree well in morphology and orientation of 
the cool spiral. In the late phase (2 Gyr), the CF in the SW 
starts to break apart in the hydro+N-body simulation. In 
the RP simulation, the CFs remain intact and the morphol- 
ogy remains close to spiral-like. We note that this break-up 
of the spiral structure can be recovered by using stronger 
damping (Fig. 6; see Sect. 3.2 for a more detailed discussion 
of the effect of damping). 


8.1.2 Size of the cold front structure 

The outwards motion of the CFs and hence growth of the 
cool spiral are important characteristics of the sloshing pro- 
cess. In their studies of the Virgo and Abell 496 clusters, 
Rlla,b‘ found that the velocity of the outward motion is 
largely independent of the subcluster, but is characteristic 
for the potential of the main cluster. This means that the 
positions of the CFs in a given cluster depend mostly on the 
time since the subcluster’s pericentre passage, i.e. the age of 
the CFs. Therefore it is important to know to what extent 
the RPA recovers this outward motion. This is the aim of 
this subsection. 

3. 1.2.1 Deriving cold front radii We study the radii 
of the CFs towards the diagonal directions in the orbital 
plane. For this purpose, we first derive radial temperature 
profiles towards NW, SW, SE and NE, where each profile is 
averaged over an azimuthal extent of ±15° (see Fig. 3 for 
profiles and Sect. 3.1.3 for their discussion). Given that CFs 
rarely form perfect circular concentric arcs around the clus- 
ter centre, this azimuthal averaging introduces a smearing 
out of the intrinsically discontinuous CFs over a finite ra- 
dial range. This occurs even for the small azimuthal range 
we use for averaging. Consequently, in the profiles, the CFs 
do not appear as a true discontinuity, but as steep slopes in 
temperature that stretch over typically 10 to 50 kpc. In each 
temperature profile, we identifv CFs as regions of tempera- 
ture slopes above 0.02 keV/ kpc. Thus, we identify an inner 
and outer edge for each CF, and the nominal CF radius is 
defined as the average radius between this inner and outer 
edge. 

The first time step at which a CF can be detected is not 
immediately after the subcluster’s pericentre passage, but 
typically 0.4 Gyr afterwards. In the hydro+N-body simula- 
tion the CFs in the NW and NE direction are established 
only at t = 1 Gyr. At even later times, in all but the NE 
direction a second CF at smaller radii is detectable. 

3. 1.2. 2 Comparing the evolution of cold front radii 

Having derived the positions of the CFs at each timestep, 
we can now proceed to analyse their outward motion. We 
do so in Fig. 2, where we plot the temporal evolution of the 
CF radii towards the diagonal directions (NE, NW, SE, and 
SW) in the orbital plane. We use error bars to indicate the 
width of each CF, i.e. from its inner to outer edge; The RP 
simulations with different settings are plotted by coloured 
lines. Here we focus on the red line, which marks our fiducial 
RP run. 

The result from the hydro+N-body simulation is the 
dashed black line. The RP results differ systematically from 


this reference: At a given timestep, the RP simulations pro- 
duce somewhat too small radii for the outermost CFs. This 
is true for all directions except NW. Thus, in general, the RP 
simulations lag behind the hydro-l-N-body one. Plotting the 
hydro+N-body result with a delay of 250 Myr (solid black 
line with error bars, and remember footnote 1) compensates 
the difference in all but the NW direction, leading to a good 
agreement to the RP simulation. The second CFs agree well 
between all runs. 

This lag in CF motion is the major systematic difference 
between two simulation methods. Consequently, we use the 
delay of 250 Myr for the full hydro+N-body simulation in 
all other comparisons regarding the fiducial merger case. 


3 A. 8 Quantitative comparison of density and temperature 
distribution 

In order to go beyond the qualitative comparison of the tem- 
perature slices in Fig. 1, we derive radial temperature and 
density profiles as described above in Paragraph 3. 1.2.1. In 
Fig. 3 we compare these profiles for different realisations for 
the fiducial merger. Again, here we concentrate on the red 
and black lines, which are for the fiducial RP simulation and 
the hydro+N-body one, respectively. 

At 0.5 GyT, the sloshing is still in the onset phase and 
we do not apply the delay discussed above. Still, already here 
the CF in the SW is ahead in the hydro+N-body simulation. 
Here, the RPA run shows a weaker impact of the subcluster 
passage on the temperature. Also the density profile in the 
northern directions are not accurately reproduced. This is 
the region the subcluster directly passes and the strongest 
differences are to be expected. 

In all later timesteps we apply the delay of 250 Myr to 
the hydro+N-body simulation as derived in Sect. 3. 1.2. 2. As 
a result, we achieve a good agreement between both meth- 
ods. Especially along the SW-NE axis, which is perpendic- 
ular to the subcluster orbit, the agreement is excellent. The 
only systematic difference between both methods is that the 
RPA over-predicts the temperature just outside the CFs, 
which we have already seen in Fig. 1. 


3.1.4 Density and temperature across the CFs 

Here we aim at comparing the evolution of the temperature 
and density at the CFs. For this purpose, we derive both 
quantities at the inner and outer edge of each CF along with 
its radius (see Paragraph 3.1.2. 1). We note that the combi- 
nation of azimuthal and radial binning introduces an uncer- 
tainty in temperature of at least ±0.5 keV at both edges. In 
Fig. 4 we plot the temperature at the inner and outer edge 
of each CF as a function of its radius. Thus, we compare 
simulations at stages when the CF spiral has reached the 
same size. We compare the results along the diagonal di- 
rections in the orbital plane. This figure is restricted to the 
outermost ring of CFs. Also here, we focus on the red and 
black lines, which are for the fiducial RP simulation and the 
hydro+N-body one, respectively. 

The temperatures at the inner edges of the CFs are re- 
produced well. The outer temperatures in the SE and SW di- 
rections differ systematically between the rigid potential and 
the hydro+N-body simulation. In these positions, the rigid 
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Figure 2. Outward motion of the cold fronts towards the diagonal directions in the ajj/-plane (one direction per panel, see label). For the 
derivation of the cold front radii, see Paragraph 3.1.2.1, The error bars indicate the width of the cold front as it appears in the azimuthally 
averaged profiles. In all but the NE panel, the evolution of the outermost and second cold front is shown. The figure compares the results 
for different realisations of the mass ratio 5 merger. The black dashed line without error bars represents the hydro+N-body run, the 
other coloured/broken lines are for the results from the rigid potential simulations with different damping settings (see legend). Clearly, 
the rigid potential simulations lag behind the hydro+N-body one. Plotting this reference with a delay of 250 Myr (and remembering 
footnote 1) corrects for this lag. 


potential simulations first underestimate the outer temper- 
ature and overestimate it at later times. The reason for this 
difference is that the rigid potential approximation leads to 
a different flow field outside the central region of the clus- 
ter. In all but the early times, this leads to a stronger com- 
pressional heating at the outer side of the CFs and thus a 
higher temperature. This systematic difference can be less- 
ened somewhat by the damping, but is not prevented com- 
pletely. 

We apply the same analysis to the density inside and 
outside each CF and present the result in Fig. 5. For the 
fiducial RP simulation, the densities at the CFs are repro- 
duced well. Thus, even at late stages, where the RPA does 
not accurately estimate the temperatures outside the CFs, 
it is still accurate for the densities. 


3.1.5 Large-scale asymmetry 

In addition to the cold spiral structure in the cluster cen- 
tre, also the large-scale distribution of the ICM density at 
least out to 500 kpc is reproduced well in the RPA. Both 
simulation methods find the characteristic asymmetry in the 
sense that profiles of, both, density and temperature on op- 
posite sides of the cluster centre alternate around each other, 
switching over at the CFs. We have illustrated this effect in 
Fig. 3 in the density panels for t = 1.5 Gyr by overplotting 
the profiles from the NE side to the SW direction and the 
NW ones to the SE direction. 

Qualitatively, the large-scale structure in temperature 


is reproduced, too. However, here the details depend on the 
damping (see Sect. 3.2). 


3.2 Impact of damping 

The damping of the inertial frame correction described in 
Sect. 2.2.3 is constructed such that it gradually switches off 
the inertial frame correction in the cluster outskirts, which 
should not be applied there. Using the RPA without this 
” re-correction” leads to an unrealistic ICM velocity field in 
the outer cluster regions and two artefacts compared to the 
hydro+N-body reference run: a cooler temperature in the 
outer northern region and the hotter temperature outside 
all CFs. Both effects can be seen in the temperature slices 
(top panel of Fig. 1) and all other comparison plots (Figs. 3, 
4 and 5). 

Using a strong damping of (Rdamp/kpc^damp/kpc) = 
(500, 300) leads to a nearly correct temperature distribution 
in the outer cluster regions as demonstrated in the temper- 
ature slice in Fig. 1. However, the temperature just outside 
the CFs is still slightly too high (Figs. 1, 3, 4). Moreover, the 
strong damping leads to slightly smaller CF radii towards 
the NW and SE direction after t = 1 Gyr, which is the di- 
rection of motion of the subcluster (Fig. 2). The CF radii 
towards the NE and SW, i.e. along the axis perpendicular 
to the orbit, are independent of damping at all times. 

In summary, the stronger the damping, the more accu- 
rate is the outer temperature distribution, but at some ex- 
pense of the accuracv in the cluster centre. Hence, we prefer 
to use a moderate damping which ensures an accurate re- 
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Figure 3. Comparison of density and temperature profiles along diagonals in the xy-plane. The columns are for different timeetcps as 
indicated in each panel. The two top rows are for the SE-NW direction, the two bottom rows for the SW-NE direction. The black solid 
lines show the hydro+N-body simulation, coloured broken lines lines the rigid potential simulations with different damping. Profiles are 
averaged over ±15° around the indicated direction. The hydro+N-body simulation is plotted with a delay of 250 Myr in all but the 0.5 
Gyr step, which accounts for the lag of rigid potential approximation described in Sect. 3. 1.2. 2. For the density profiles at t = 1.5 Gyr, 
we overplot the profile of the positive direction as thin lines in the negative direction in order to demonstrate the asymmetry. 


production of the CF radii, temperatures inside them and 
densities at, both, the inside and outside. 

The case (800, 100) is very similar to our fiducial setting 
(800, 300) demonstrating that the results are not sensitive 
to the choice of the fall-off length scale for the damping, 

^damp* 


3.3 Subcluster orbit and mass evolution 

In addition to the ICM properties, we also compare the sub- 
cluster orbit w.r.t. the cluster centre from both methods in 
Fig. 7. The test particle orbit used in the RP simulations ap- 
proximates the true trajectory of the subcluster very well. 
We have marked the subcluster position in steps of 250 Myr 
along both orbits, demonstrating that the test particle ap- 
proximation also recovers the motion of the test particle 
along its orbit up to a few 100 Myr after pericentre pas- 
sage. After that, dynamical friction causes the subcluster 
to decelerate significantly, such that its new apocentre dis- 
tance is only about 1.8 Mpc, and the apocentre is reached 
already 1.25 Gyr after pericentre passage. The test particle 
method does not capture this effect, and predicts a com- 
parable cluster-centric distance already after about 0.6 Gyr 
after pericentre passage. Taylor &; Babul (2001) proposed 


an algorithm to incorporate dynamical friction in a simpli- 
fied form. However, the gas sloshing we are interested in is 
triggered by the pericentre passage and thus is unaffected 
by this difference in subcluster position at later epochs. 

In Fig. 8 we demonstrate the evolution of the subcluster 
by plotting its mass within its scale radius. Clearly, during 
pericentre passage and up to 0.2 Gyr afterwards the sub- 
cluster suffers tidal compression. It starts loosing mass sig- 
nificantly only 0.5 Gyr after pericentre passage. This is well 
after triggering the gas sloshing and thus has no significant 
effect on the further evolution of the ICM in the main clus- 
ter. We will discuss the differences in the evolution of the 
main cluster potential in Sect. 4.1. 

3.4 Mergers with other mass ratios 

8.4 1 Mass ratio 20 

We have run the same comparison for a minor merger with 
a larger mass ratio of 20. Here, the subcluster has a smaller 
scale radius of only 220 kpc, hence the damping setting 
(i^damp/ kpc, Tdamp/ kpc) = (500, 300) corresponds to mod- 
erate damping here. In Fig. 9 we compare temperature slices, 
in Fig. 10 the evolution of the CF radii, and the tempera- 
tures inside and outside each CF in Fig. 11. Even for this 
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Figure 4 Temperature inside and outside of the cold front as a function of cold front radius. We show results for the four diagonal 
directions in the a^-plane. Different line colours/styles code different realisations of the mass ratio 5 merger, see legend. In each panel, 
the upper and lower sets of lines denote the temperatures just outside and inside the cold fronts, respectively. For comparison, we plot 
the initial temperature profile as the dotted line. 



Figure 5. Same as Fig. 4, 


50 100 150 200 250 50 100 150 200 250 

r CF f kpe 

but for density inside and outside of the cold fronts (upper and lower set of lines, respectively). 


case with a much smaller subcluster, the sloshing in the 
RPA runs lags behind the hydro+N-body one by 200 Myr. 
Hence, ve come to the same results and conclusions as in 
the fiducial case. 


3.4-2 Mass ratio 2 

A merger with a mass ratio of 2 is no minor merger any- 
more, and we perform this simulation with the purpose of 


exploring the limits of the RPA. Here, the subcluster has 
a scale radius of 416 kpc and we use a damping setting of 
(-Rdamp/kpc, Ldamp/kpc) = (950,400). In Fig. 12 we com- 
pare temperature slices, in Fig. 13 profiles in the orbital 
plane along the diagonal directions at different timesteps. 
For this case, we find a delay of 250 Myr between RPA and 
hydro+N-body ones. 

In both methods, the central ICM starts sloshing in the 
cluster centre, and a clear primary CF is formed in the SW. 
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Figure 12. Snapshots of the temperature in the orbital plane for the merger with mass ratio 2. We show the hydro+N-body results (upper 
row) with a delay of 250 Myr compared to the rigid potential run (bottom row). This hydro-|-N-body simulation has a low resolution of 
10 kpc, which is at least partially responsible for the loss of the cool dense centre. The major cold front in the SW is still reproduced 
well, but T-ie general distortion of the cluster by this major merger is not captured anymore by the rigid potential approximation. 



i 1 1 — i 1 1 — i 1 1 1 t ,.ii 1 1 i . — a. i i 
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Figure 13. Comparison of density and temperature profiles along the diagonals in the xy-plane. Same as Fig. 3 but for a mass ratio of 
2 between the clusters, which is already a major merger. The hydro+N-body simulation is plotted with a delay of 250 Myr. 
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Figure 6. Temperature slices for the fiducial run at t = 1.5 Gyr, 
for rigid potential simulations without damping (upper panel) 
and with strong damping (lower panel). Compare to Fig. 1 for 
moderate damping and hydro+N-body case. 



x/kpc 

Figure 7. Comparison of the subcluster orbit in the hydro+N- 
body simulation (black solid line) and the rigid potential approx- 
imation (red dashed line) for the merger with mass ratio 5. We 
mark the subcluster positions in 250 Myr intervals. While the 
direction of the orbit is reproduced well in the rigid potential ap- 
proximation, the neglect of the dynamical friction leads to signifi- 
cant over-estimation of the velocity after 0.5 Gyr after pericentre 
passage. 


However, in the hydro+N-body code the cool central gas 
moves completely out of the cluster centre, forming a CF 
towards the NE of the centre only at very late stages. In 
the RPA, the cool gas core is never completely displaced 
from the potential minimum. The comparison of profiles in 
Fig. 13 reveals that the RPA still achieves a good agreement 
for the primary CF in the SW and SE, despite the nearly 
equal masses of both clusters. 

Thus, even for this major merger, the RPA recovers the 
evolution of the major CF. However, it comes to its limit 
regarding the morphology, the other CFs, and the large-scale 
distortion of the main cluster. 



Figure 8. Evolution of subcluster mass within the subcluster 
scale radius. During pericentre passage and up to 0.2 Gyr after- 
wards the subcluster suffers tidal compression. It starts loosing 
mass significantly only 0.5 Gyr after pericentre passage. 
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Figure 9. Snapshots of the temperature in the orbital plane for 
the merger with mass ratio 20. We show the hydro-fN-body re- 
sults (upper row) with a delay of 200 Myr compared to the rigid 
potential run (bottom row). Both methods give very similar re- 
sults. 

4 DISCUSSION AND SUMMARY 

4.1 Origin of the delayed evolution in the RPA 

We have shown that the rigid potential approximation can 
reproduce the characteristics of gas sloshing well except for 
two artefacts, the higher temperatures at the outside of the 
CFs and the temporal lag in evolution. We have traced back 
the former effect to the necessarily unrealistic gas motions 
in the outer cluster region in this approximation. In order to 
investigate the origin of the latter, we compare the evolution 
of radial potential profiles in the fiducial case in Fig. 14. 
We show profiles towards the N, S, W and E, where we 
plot one direction per panel. Different timesteps are colour- 
coded (coded by line style in print version). The results from 
hydro+N-body and from the RPA are shown by solid lines 
and dashed lines, respectively (thick and thin lines in print 
version). 
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Figure 14. Evolution of radial potential profiles in the main cluster towards N, W, S, E. The profiles are averaged azimuthally over ±5°. 
Hydro+N-body simulations are shown by solid lines (thick lines in print version), rigid potential simulations by dashed lines (thin lines 
in print version). The colours (line styles in print version) code different timesteps, see legend. In both methods, the subcluster passage 
leads to a temporary deepening of the central potential. Additionally, the hydro+N-body method captures the tidal compression of the 
main cluster centre, leading to a stronger deepening and steepening of the potential. Both effects wear off with time and 1 Gyr after 
later the initial potential is nearly recovered. However, already after about 0.25 Gyr, the slopes of the potential are very similar in both 
methods. 




Figure 10. Same as Fig. 2, but for the merger with a mass ratio Figure 11. Same as Fig. 4, but for the merger with a mass ratio 
of 20. The result from the hydro+N-body simulation is plotted of 20. 
with a delay of 200 Myr. 


Generally, the central potential deepens and steepens 
during the passage of the subcluster. The direct overlap with 
the subciuster potential can also lead to a temporary local 
flattening of the potential, e.g. at t =* 0.1 Gyr in the E. 
In the RPA, the effect is solely due to the overlap of both 
potentials. The hydro+N-body code additionally captures 
the tidal compression of the cluster centre, which leads to 
the extra deepening and steepening during the pericentre 


passage evident from Fig. 14. This leads to a temporarily 
different evolution of the central potentials between both 
methods. At 1 Gyr after pericentre passage, in both methods 
the potential is nearly back to its initial state, where the 
closest ” recovery” is achieved in the S. This explains why 
the major, southern CF evolves so similar in both methods. 

In order to access the timescales of the potential evolu- 
tion more directly, we study the evolution of the radial gravi- 
tational acceleration in different positions. Along each direc- 
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Figure 15. Evolution of radial gravitational acceleration at different radii towards N, W, S, E. Comparison between hydro+N-body 
simulations (solid lines) and rigid potential simulations (dashed lines). The colours code different radii, see labels. The effect of the 
subcliistei passage is clearly seen in both methods. The effect is stronger in the hydro+N-body method because it captures the tidal 
compression of the cluster core. After about 0.25 to 0.5 Gyr, the acceleration is again similar in both methods. Only in the inner 100 
kpc, the acceleration is slightly lower in the hydro+N-body method. Thus, after a slightly different onset, the sloshing proceeds in a very 
similar way afterwards. 


tion, we calculate the gravitational acceleration at four radii 
and plot their temporal evolution in Fig. 15. Again, there is 
one panel per direction, and the radii are colour-coded. Solid 
lines are for the hydro+N-body code, dashed for the RPA. 
Here we see more clearly the modification of the gravita- 
tional field described above. Prior to pericentre passage, the 
accelerations in both methods agree well. During pericentre 
passage, the tidal compression leads to stronger accelerat- 
ions in the hydro+N-body method. After about 0.25 to 0.5 
Gyr, the accelerations have reached the final levels at which 
they remain until the next pericentre passage. Outside 200 
kpc, the initial acceleration is nearly recovered, whereas in- 
side 100 kpc the final acceleration is slightly lower than the 
original one. The similarity of the accelerations, i.e. the po- 
tential slopes, after the onset of sloshing accounts for the 
similar evolution of the CFs in both methods. The period of 
different potential evolution at and shortly after the pericen- 
tre passage is responsible for the different evolution during 
this phate and sends the hydro+N-body method ahead, be- 
cause in a steeper potential the sloshing oscillations tend to 
be faster. 


4.2 Choice of damping parameters 

In Sect. 2.2.3 we suggested that the damping radius, ftdamp, 
outside 'vhich the inertial frame correction will be switched 


off, should be comparable to the diameter of the subclus- 
ter, i.e. twice the subcluster scale radius. We investigated 
the impact of the damping parameters (Sect. 3.2) for all 
cases studied here and found that this choice leads indeed 
to the best possible reproduction of all features. Moreover, 
we find that the evolution of CF radii is the same with this 
mild damping and in the undamped case, but CF radii are 
somewhat lower if the damping is chosen too strong. This 
difference in behaviour can be used to confirm the choice of 
f^damp for clusters where no hydro+Nbody simulation is yet 
available. As long as the CF radii in the damped simulation 
agree with the ones in the undamped case, the choice of 
Rd&mp is reasonable. Regarding the damping scale length, 
f'damp, we found that our results are not sensitive to this 
parameter, and we recommend using Ldamp « 0.5f?damp, as 
we have done here. We also followed this strategy for choos- 
ing iidamp and Ldamp in our simulations of Virgo and A496 
(Rlla,b). 


4.3 Summary - reliability of the rigid potential 
approximation 

We investigated the reliability of the rigid potential approxi- 
mation described in Sect. 2.2 for simulations of minor merger 
induced gas sloshing. We use the hydro+N-body simulations 
of Z10 of the same scenario as the reference. Those capture 
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the full evolution of the ICM and DM components including 
dynamical friction, tidal compression and tidal stripping. In 
contrast, the rigid potential approximation treats the poten- 
tials of the individual clusters as static and models only their 
relative motion. This simplification makes the rigid potential 
simulations faster by about a factor of 5 for the resolution 
used here, and more for higher resolution due to the im- 
perfect scaling of the Poisson solver. This speed-up is very 
useful in several circumstances. For example, constraining 
the merger history of a given cluster by reproducing the ob- 
served sloshing signatures usually requires a large set of sim- 
ulations (see Rlla for Virgo cluster and Rllb for Abell 496). 
Also investigations of the impact of more time-consuming 
physics like viscosity (Z10) and thermal conduction benefit 
from a fast method for the basic process. 

While we expect and find temporal differences in the 
evolution of the gravitational potentials, regarding the gas 
sloshing, our main interest lies in the evolution of the ICM. 
Hence, we have simulated three representative merger sce- 
narios using both methods and compared them in detail. We 
have shown that, except for two (correctable) artefacts, the 
rigid potential approximation reproduces the results of the 
hydro+N-body runs very well: 

(i) The rigid potential approximation reproduces the typ- 
ical sloshing cold fronts and central cold spiral structure in, 
both, morphology and orientation (Figs. 1, 9, 12). 

(ii) The minor merger induces a characteristic large-scale 
asymmetry in the main cluster’s ICM beyond the central 
cool spiral structure, which is reproduced at least qualita- 
tively in, both, density and temperature (Fig. 3). 

(iii) The outward motion of the cold fronts with time is 
delayed in the rigid potential approximation by typically 
200 Myr compared to the hydro+N-body one. This origi- 
nates from a temporarily different potential shape during 
pericentre passage because the rigid potential approxima- 
tion does not capture tidal compression. The significantly 
different evolution is, however, short-lived, and the evolution 
proceeds very similar afterwards. Therefore, if this delay is 
corrected for, also the outward motion of the cold fronts is 
reproduced very well by the rigid potential approximation 
(Figs. 2, 10). 

(iv) The ICM density on both sides of the cold fronts is 
reproduced well, and so is the temperature inside the cold 
fronts (Figs. 4, 5, 11). The temperature at the outside of the 
cold fronts tends to be slightly too high compared to the ref- 
erence simulations due to a necessarily unrealistic velocity 
field in the cluster outskirts in the rigid potential approxi- 
mation (Sect. 3.2). This effect is strongest if only the basic 
inertial frame correction is used (Sect. 2.2.2), and is milder 
for the improved version, where the inertial frame correction 
is only applied to the central part of the galaxy cluster and 
gradually dampened towards large radii (Sect. 2.2.3). Too 
strong dampening, however, causes slightly too small cold 
front radii towards the directions aligned with the subclus- 
ter orbit. Our tests recommend dampening outside cluster- 
centric r:dii of about twice the subcluster scale radius, with 
a characteristic fall-off scale length comparable to the sub- 
cluster scale radius. 

Thus, the rigid potential approximation method can be em- 
ployed e.g. in order to disentangle the merger history for an 
observed cluster (Rlla,b). The agreements in items (i) and 


(ii) guarantee that the orientation of the subcluster orbit 
can be inferred correctly. Item (iii) ensures that the age of 
the cold front can be estimated reasonably, although here the 
uncorrected rigid potential approximation will over-estimate 
the age by about 200 Myr. If the age needs to be determined 
more accurately, especially during the onset of sloshing or 
the early evolution of the cold fronts, hydro+N-body sim- 
ulations are required. The mass of the subcluster can in 
principle be decoded from the contrast of density and tem- 
perature across the cold fronts. However, this attempt is in- 
trinsically difficult, because usually the cold fronts are found 
within the cool cores of their host clusters, where the gen- 
eral gradient of, both, the density and temperature is in the 
same direction as the cold front discontinuity itself. Given 
that azimuthal and radial averaging in deriving radial pro- 
files smears out the cold fronts over a certain radial extent, 
the contrast of all quantities between the inner and outer 
edge of a given cold front includes the general variation of 
the quantity across in this radial range. The contributions 
of the general profile and the cold front discontinuity to this 
contrast are hard to separate, no matter which simulation 
method has been used. 


ACKNOWLEDGMENTS 

ER is supported by the Priority Programme ’’Witnesses of 
Cosmic History” of the DFG (German Research Founda- 
tion) and the supercomputing grants NIC 3711 and 4368 
at the John-Neumann Ins tit ut at the Forschungszentrum 
Jiilich. JAZ is supported under the NASA postdoctoral pro- 
gram. We thank Marcus Briiggen for helpful discussions, 
and the referee Max RufFert for his clarifying comments. 
The results presented were produced using the FLASH code, 
a product of the DOE ASC/Alliances-funded Center for 
Astrophysical Thermonuclear Flashes at the University of 
Chicago. 


REFERENCES 

Ascasibar Y., Markevitch M., 2006, ApJ, 650, 102 

Bourdin H., Mazzotta P., 2008, A&A, 479, 307 

Churazov E., Forman W. R., Jones C., Bohringer H., 2003, 
ApJ, 590, 225 

Clarke T. E., Blanton E. L., Sarazin C. L., 2004, ApJ, 616, 
178 

Dubey A., Antypas K., Ganapathy M. K., Reid L. B., Riley 
K., Sheeler D., Siegel A., Weide K., 2009, Parallel Com- 
puting, 35, 512 

Dupke R., White III R. E., Bregman J. N., 2007, ApJ, 671, 
181 

Fabian A. C., Sanders J. S., Taylor G. B., Allen S. W., 
2005, MNRAS, 360, L20 

Ghizzardi S., Rossetti M., Molendi S., 2010, A&A, 516, 
A32 

Hemquist L., 1990, ApJ, 356, 359 

Johnson R. E., ZuHone J. A., Jones C., Forman W., Marke- 
vitch M., 2011, eprint arXiv: 1106.3489 

Markevitch M., Gonzalez A. IL, David L., Vikhlinin A., 
Murray S., Forman W. R., Jones C., Tucker W., 2002, 
ApJ, 567, L27 



14 E. Roediger & J. ZuHone 


Markevitch M., Ponman T. J., Nulsen P. E. J., Bautz 
M. W., Burke D- J., David L. P., Davis D., Donnelly R. H., 
l<brman W. R., Jones C., Kaastra J., Kellogg E., Kim D., 
Kolodziejczak J., Mazzotta P., Pagliaro A., Patel S., Van 
Speybroeck L., Vikhlinin A., Vrtilek J., Wise M., Zhao P., 
2000, ApJ, 541, 542 

Markevitch M., Vikhlinin A., 2007, Physics Reports, 443, 
1 

Markevitch M., Vikhlinin A., Mazzotta P., 2001, ApJ, 562, 
L153 

Mazzotta P., Edge A. C., Markevitch M., 2003, ApJ, 596, 
190 

Mazzotta P., Giacintucci S., 2008, ApJ, 675, L9 
Mazzotta P., Markevitch M., Vikhlinin A., Forman W. R., 
David L. R, VanSpeybroeck L., 2001, ApJ, 555, 205 
Million E. T., ALlen S. W., 2009, MNltAS, 399, 1307 
Owers M. S., Nulsen P. E. J., Couch W. J., Markevitch M., 

2009, ApJ, 704, 1349 

Roediger E., Briiggen M., Simionescu A., Bohringer H., 
Churazov E., Forman W. R., 2011, MNRAS, 413, 2057 
Roediger E., Lovisari L., Dupke R. A., Ghizzardi S., 
Briiggen M., 2011, MNRAS, submitted 
Sanders J. S., Fabian A. C., 2006, MNRAS, 371, 1483 
Sanders J. S., Fabian A. C., Dunn R. J. H., 2005, MNRAS, 
360, 133 

Sanders J. S., Fabian A. C., Taylor G. B., 2009, MNRAS, 
396, 1449 

Simionescu A., Werner N., Forman W. R., Miller E. D., 
Takei Y., Bohringer H., Churazov E., Nulsen R E. J., 

2010, MNRAS, 405, 91 

Taylor J. E., Babul A., 2001, ApJ, 559, 716 
Vikhlinin A., Markevitch M., Murray S. S., 2001, ApJ, 551, 
160 

ZuHone J. A., Markevitch M., Johnson R. E., 2010, ApJ, 
717, 908 

ZuHone J. A., Markevitch M., Lee D., 2011, eprint 
arXiv: 1108.4427 

This paper has been typeset from a LfeX/ I£I]eX file prepared 
by the author. 



